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p Abstract 

A new formula for the heavy quark-antiquark spin dependent potential is given 
by using the techniques developed in the heavy quark effective theory. The leading 

x 

5_i ■ logarithmic quark mass terms emerging from the loop contributions are explicitly 

extracted and summed up. There is no renormalization scale ambiguity in this new 
formula. The spin-dependent potential in the new formula is expressed in terms of 
three independent color-electric and color-magnetic field correlation functions, and it 
includes both the Eichten-Feinberg's formula and the one-loop QCD result as special 
cases. 
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The study of the structure of the spin- dependent (SD) interaction potential between 
heavy quark and antiquark from QCD is one of the interesting problems in heavy quark 
physics. So far there are mainly two kinds of approaches in the literature. The first kind 
of approach starts from the static limit (infinitely heavy quark limit) and makes relativistic 
corrections via the 1/m expansion, where m stands for the heavy-quark mass. The for- 
mula for the SD potential in this approach was first given by Eichten and Feinberg [[jj] in 
which the potential is expressed in terms of certain correlation functions of color-electric 
and color-magnetic fields weighted by the Wilson loop factor. Later Gromes [f| derived 
an important relation between the spin-independent (SI) and the SD potentials from the 
Lorentz invariance of the total potential and the correlation functions given in Ref. JIJ,and 
it supported the intuitive color-electric flux tube picture of color-electric confinement sug- 
gested by Buchmiiller . The second kind of approach is to calculate the SD potential from 
perturbative QCD up to one-loop level and put in the nonperturbative part of the poten- 
tial by hand @||. In this approach, certain terms containing lnm emerge from the loop 
contributions. Furthermore, in the case of unequal masses, new structure of the order-l/m 2 
spin-orbit coupling containing In m arises in this approach, which is not included in the first 
kind of approach It seems that there exists a discrepancy between these two kinds of 

approaches ||. To understand the essence of this discrepancy and to have a deeper under- 
standing of the In m-dependence in the SD potential has become an important theoretical 
problem in heavy quark physics for a quite long time. In a recent paper [[Fj, two of us (Chen 
and Kuang) derived some general relations between the SI and SD potentials in the spirit 
of the 1/m expansion by using the technique of the reparameterization invariance 



developed in the heavy quark effective theory (HQET), which include the Gromes relation 
and some new useful relations, and lead to the conclusion that the general structure of the 
SD potential is the same as that obtained in the second approach. However, such a simple 
symmetry argument does not concern the problem of the In m-dependence of the potential. 

In this paper, we adopt the conventional formulation in the HQET to construct an 
effective QCD Lagrangian including both a heavy quark field and a heavy antiquark field, 
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with which we study the SD interaction potential. The conventional HQET has proved to 
be very powerful in studying the heavy-light quark systems pT^p^ pl HT^jr^ri^l and most of 



the useful techniques developed in it can be applied to the present theory We emphasize 
here that we are working at the quark level to study the interaction potential between heavy 
quarks rather than studying the bound state quarkonia. The reason causing the discrepancy 
between the two approaches can be easily understood in the effective Lagrangian formalism. 
When an effective Lagrangian is constructed by expanding the full theory in terms of 1/m, 
the ultraviolet-behavior of the theory is changed. Therefore, when the order of the 1/m 
expansion and the loop integration are exchanged, which is involved in passing from the full 
theory to the effective theory, differences such as logarithmic quark mass terms can emerge. 
In order to reproduce the results of the full theory, one must match the full theory and 
the effective one. With this matching condition these logarithmic quark mass dependent 
terms can be explicitly extracted in the coefficients of higher dimensional operators and can 
be summed up by using the renormalization group equation (RGE). In the original work 
of Eichten and Feinberg |I[ the authors expanded the full heavy fermion propagator in an 
"external" gluon field A 11 in terms of 1/m. It is easy to see that each term in their expression 
corresponds to insertions of operators in the effective Lagrangian which contribute to the SD 
potential up to order 1/m 2 , but with only tree level coefficients. Therefore, the logarithmic 
terms are not accounted for. This is the essential reason that there is a discrepancy of 
the quark mass dependence between the Eichten-Feinberg-Gromes (EFG) formula and the 
one-loop calculation. With this insight, we can improve the EFG formula by summing 
the logarithmic quark mass terms and establish a consistent picture reconciling these two 
kinds of approaches in the framework of the effective Lagrangian. To this end, we need not 
only to renormalize the effective Lagrangian up to order of 1/m 2 , but also to consider the 
mixing of the nonlocal operators with local four fermion operators when we use the HQET 
to calculate the heavy quark-antiquark Green's functions. Then we can follow the methods 
developed in Ref. P~5|JT[] and obtain a new formula for the spin dependent potential by using 
the renormalized effective Lagrangian. The above results of the two approaches are just two 



special approximations of our new formula. Moreover, our new formula is independent of 
the renormalization scale parameter /i, so that it does not suffer from the scale ambiguity 
as the second approach does. 

Let us first construct the renormalized effective Lagrangian for one heavy quark field up 
to order of 1/m 2 . We start from the full QCD Lagrangian. As in the conventional HQET, 
we define the heavy quark field h v+ (x) and the heavy antiquark field h' v _{x) related to the 
original field ip(x) as 



h v+ (x) = P + e imv ^(x), h' v _(x) ee P_e- irmMS iP(x), (1) 

1 ±4 

where v is the velocity of the heavy quark, and P± ee — - — . Integrating out the quantum 



fluctuation of the quark field, we obtain the heavy quark effective Lagrangian C c [|l"2| , |10l . 
After expanding order by order in powers of 1/m, the first three terms are 

C = c h v+ (x)iD ■ vh v+ (x), (2) 

Ci = c$cih v+ (x)— h v+ [x) - c c 2 h v+ {x) — h v+ {x) + c c 3 g s h v+ {x) — h v+ [x), (3) 

2m 2m Am 

^ = c 01 ^ l h v+ (x)(c 4 v v Di>G IJV + ic b a^v°D^G va )h v+ {x) 



'4m 2 
Am z 



cSDf-ic 7 {D-v) 2 -c^G^ 



iD ■ vh v+ (x), (4) 



where iD^ = id^ — gA a ^T a . The last term in £2 has no contribution in order 1/m 2 due to 
the equation of motion. At tree level the q's are all unity. Note that after the expansion 
the high energy behavior is different from that in the full theory. So the operators in (@)-(fD 
need to be renormalized and their coefficients can be determined by matching to the full 
theory. Here ^Jcq corresponds to the wavefunction renormalization constant, and c\ — C3 



have been calculated in Refs. JL31 and |14], using the RG summation, Ref. [T^ gives 



_9_ 
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ci{fjt,m) = l, c 2 (/i,m) = 3 -2, c 3 (/i,m) = , (5) 

\a s (m)J \a s {m)J 

where q ee Cj(/i, m). The coefficients C4 — cq can be determined by the reparameteriza- 
tion invariance of C c which leads to the reparameterization invariance of the renormalized 
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Lagrangian. This is shown as follows: Consider the infinitesimal velocity transformation 



v — > v + Av pi. The infinitesimal transformation of 5C can be written as 



S£ = 5% + ^-STx, (6) 
2m 



where 



S% = c (l ~ c 1 )h v+ (x)iD ■ Avh v+ (x), (7) 



and 



5Ti = co(l — 2c 2 + cq) h v+ (x)iD ■ AviD ■ vh v+ (x) 

+ c (c 3 - c 2 + ~c 4 - ~c 5 ) h v+ (x)Yv u igG '^^h v+ (x) (8) 
+ c (l - c 2 - c 3 + ic 4 + |c 5 ) h v+ (x)^Yv v igG^h v+ (x). 

Different terms in (0)-(^D are of different Lorentz structures, and therefore, they should 
vanish separately. 5Tq = leads to C\ — 1, i.e., the kinetic energy term is not renormalized 
0, and 5Ti = gives the two relations 

c 4 = c 6 = 2c 2 - 1, c 5 = 2c 3 - 1. (9) 

Using Eq.(|) we obtain 

_ JL _ JL 

c 4 (^m) = 6f^4) " -5, c 5 (/i,m) = 2f^4') " - 1. (10) 
\a s (m) J \a s {m) J 

The effective Lagrangian for antiquark field can be obtained by simply replacing v by 
— v and h +v (x) by h'_ v {x) in the above effective Lagrangian. 

To study the heavy quark-antiquark interaction potential we are going to evaluate the 
heavy quark-antiquark four point Green's function to order 1/m 2 . In doing this, we need to 
calculate all the possible order 1/m operator insertions. Then, additional divergences will 
appear from double insertions of these operators, i.e., these bilocal operators will mix with 
certain local four-fermion operators. Let us consider a general unequal mass case, e.g. the 
heavy quark Qi with mass nil and the antiquark Q 2 with mass m 2 , as in the case of the 
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cb. To order 1/m 2 , there are only two local dimension 6 color singlet four fermion operators 
Oi(x) and 2 (x): 

9 2 - 

Oi(x) = - — - — h v+i (x)a^h v+j (x)h' v - j (x)a fMl/ h' v _ i (x), (11) 

g 2 - 

2 (x) = - — - — h v+i (x)a flu h v+ i(x)h' v - j (x)a l , l/ h' v _ j (x), (12) 

where i,j = 1,2, • • ■ ,N C are color indices. Suppose m% > m 2 . The heavy quark antiquark 
effective Lagrangian up to order l/mf, l/m^, and 1/m 2 , can be constructed in two steps 
as follows. Starting from the Lagrangian in the full theory, we first treat Q\ as a heavy 
quark, and obtain 

We do not need to add new operators in (|13|) because there are no divergent terms of the 
form of Oi(x) and 2 (x). Next, we treat Q 2 as a heavy antiquark, and obtain 

£" = £-Qieff + C-Q-zeff + rfl(/i)Oi(/i) + d 2 (/J,)0 2 (/J,), (14) 

where the last two terms are the two necessary dimension-6 operators with unknown coeffi- 
cients di(fi) and d 2 (n), respectively. Now we determine di(fi) and d 2 ([i) by using the RGE. 
It is easy to see that only the magnetic operator insertion in each fermion line will mix with 
0\ and 2 due to the Lorentz structure. Let us denote the resulting contribution by 

2 

O (x) = — ^ / d 4 yT* \K+{x)G^h v+ {x)h' v _{y)G^h! B _{y)] , (15) 

and its coefficient by do(/j,). Here T* means time ordering. The coefficients d(fi) = 
(do(fi), di(fi), d 2 (n)) satisfy the renormalization group equation 

^d(fi) + d( f i) 1 = 0, (16) 

where 7 is the anomalous dimension matrix. A straightforward one-loop calculation using 
the HQET technique gives 
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4vr 2 



-N r — 



1 \ 



V 



• (17) 


means that local operators are not able 



Here 7 00 = 27 mas . The results 710 = 720 
to mix with bilocal operators. However, bilocal operators can mix with local operators, 
and there are two box, two cross, and a "fish" diagram contributing to 701 and 702. That 
7n = 712 = 721 = 722 = can be understood as follows: The anomalous dimensions are 
gauge independent. If we take an axial gauge where v ■ A = gluonic interactions decouple 
from the fermion field in the zeroth order effective Lagrangian (0). The initial condition for 
the RGE is determined by matching the effective theory to the full theory at \i — m 2l i.e. 



d{m 2 ) = ( c 3 (m 2 , mi), 0,0 ) . 



;is) 



With this, the solution of the RGE ([16]) are 

( a 2 s{») 
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do(n) = c 3 (/i, m 2 )c 3 (/i, mi) 
1 



d 2 (/i) 



c 3 (m 2 ,mi)[l - c 3 (/i,m 2 )] 
1 



a s (m 1 )a s (m 2 ) J 

__o_ 

1 I a,(m 2 )\ 25 



a s K) / 



18 
" 25 



(19) 



N, 



■di(n). 



Our effective Lagrangian ( 14[) is thus completely determined. 

Now we apply it to calculate spin-dependent force. We shall take v = (1,0,0,0) and 
denote h v+ (x), h' v _{x) as h(x), h'{x) for short. Similar to Ref. we introduce a gauge 
invariant four-point Green's function 

/ = (0\T*[h'(y 2 )f B P(y 2 ,y 1 )h(y 1 )] [^(arOr^P^!, x 2 )h'{x 2 )\ |0>, (20) 

where P(x, y) = Pexp ig J* dz^A^^z) is the path-ordered exponential |IB| , |I7| . As is argued 
in Refs. |15|,[TJ, in the limit that the time interval T = (y® + y 2 )/2 — (x? + x 2 )/2 — > 00, with 
x\ — x\ and y\ — y\ fixed, the limit of / is 



I 6 AB 6(r x - r y ) exp[-Te(r) 



(21) 



where r x = x\ — x 2 , r y — yi — y 2 , r = \ f x \, and e(r) is just the static energy between the 
quark and antiquark separated by the spatial distance r. Here the appropriate ordering of 
the limits is that first m — > oo and then T — > oo, so that the motion of the quark and the 
antiquark can be treated perturbatively [0,0]. 

Taking all the operators with dimension higher than 4 in the Lagrangian as the pertur- 
bative part, / can be calculated by using standard perturbation theory. In the calculation, 
the zeroth order full fermion propagator So(x,y,A) in the external gluon field is used 
and it is PJ| 



S (x, y, A) = -id {x - y )P(x , y Q )5(x - y). 



Next we define the symbol (• • •) and / as 



(22) 



exp(iS YM (A)), (23) 



I = Tr{P + Y P_r/)5(f ! - yi)5(f 2 - y 2 ). 

To order 1/m 2 , l/(mim 2 ) and 1/m-i, / can be expressed as 
rT/2 1 

dz — (D (■k 1 ,z) - Cz{yi,mx)g a {yL)a 1 ■ B(-x h z)) + (1 <-> 2) 
t/2 rrii 



I = 

9s(») f T ' 2 



4mf Jt/2 

I rT/2 /-T/2 



dz[(c 4 (/i,mi)^- - c 5 (//,mi)iei ii; cr^)(P l (xi,z)D ; '(xi,2;)) + (1^2) 



/ dz dz'8(z' - z)[((D 2 - c^ mi )g s {^)a\B % ){^z) 

nil J -T/2 J -T/2 

1 rT/2 rT/2 

(D 2 - c 3 (fx,mi)g s ((x)a\B l )( X i,z')) + (1 <-> 2)] / cfe / dz' 

mini2 Jt/2 J-t/2 

((D 2 - c 3 (/i,mi)^(/x)ai J B i )(xi,z)(D 2 - g s (u)c^ ) m 2 )a 3 2 B j ){x 2l z')) 



+ ^^Td^)a\a\5\-Ki-x 2 ), 



2niim 2 



where 



d(fi) = di(u) + 



c^) 



(24) 



(25) 



(26) 



Similar to the derivation given in Ref. [J]], we obtain the spin-dependent potential 

1\ dV (r) 



V(r) = V (r) + ( — ^ + — | I -L 



where 



Si S2 

2 2 
mf W2 



c + (fi,m 1 ,m 2 ) 



dr 



+ c + (fi,m 1 ,m 2 ) 



dr 



[ S 1 + S 2 \ ,l dV 2 {fx,r) 

+ ■Lc + (fi,m 1 ,m 2 ) 

V m\vn 2 I r dr 



1 



+ 



+ ^ 



1 (Si ■ r)(S 2 ■ r) - -Si ■ S 2 

m\m 2 r 2 
1 1 



c 3 (/i, mi)c 3 (/i, m 2 ) V^, r) 



Si ■ S 2 (c 3 (/i, mi)c 3 (/i, m 2 )\/ 4 (/i, r) - 6N c g 2 s (n)d(n)5(r) 



3 m,\m 2 

(Si S 2 \ ,l rf[%,r) + yi(/z,r)] 

+ — 2 -Lc-O^m^ma) — 1 

\ mf mn / r dr 



+ 



Si-S 2 \ 

m\m 2 I 



ldV 2 (fx,r) 
•Lc_(/i, m 1 ,m 2 )~- 



r dr 



(27) 



c + (ii,m 1 ,m 2 ) = ^[c 3 (/i,mi) + c 3 (/i,m 2 )], c_(//,mi,m 2 ) = ^[c 3 (/i,mi) - c 3 (/i,m 2 )], (28) 



and 
Vo(r 



Urn ' n<1> 



T^oo T 

IdV^r) 



r-T/2 pT/2 I z ' _ z 

r k -— — = lim e ijk / dz dz' —— 

r dr t^oo J J-t/2 J-t/2 \ T 



gl{n)l2{B\^z)E^^z'))l{l), 



1 dVo(u r) r T / 2 r T / 2 I z'\ 

r k - 2 7' ' = lim e ijk / dz dz' ^ ^(/i)/2(^(x 2 , z)^(x lf z')>/<1>, 

r Ctr T^oo J-T/2 W / 



(29) 
(30) 

(31) 



5 ij S ij r T l 2 r T / 2 a 2 (u) 

Wj - -)V^r) + °—Vi(jM,r)] ee lim / dz / dJ^(B>(x 1 ,z)B>(x 2 ,J))/(l). (32) 

3 3 T-»ooJ_r/ 2 J-T/2 J 

In Ref. 0, Gromes derived a relation J~l%( r ) + ^ / i(/ i ? r ) — V 2 (fj,, r )] = 0. The relation 
c§{n,m) = 2c 3 (/i, m) — 1 in @ obtained from reparameterization invariance ensures that 
the general relations derived in Ref. are all satisfied. It also shows that those relations 
are consistent with each other. Using those relation, the SD potential can be simplified as 
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V^(r) = Vb(r) + ( — H 5 ) -L ( c+(//, mi, m 2 ) 



v mf m2 J \ dr 2 dr 

( Si + S 2 \ ,ldV 2 {n,r) 
+ •Lc + (^,rai,ra 2 ) 



m\m 2 / r dr 

x (Si-r)(S 2 -r)-^Si-S 2 r 2 

5 — c 3 (/i, mi)c 3 (/i, m 2 )y 3 (/i, r) 

mim 2 

1 1 



(33) 



3 m\m 2 

+ 



Si ■ S 2 (c 3 (/i,mi)c 3 (/i,m 2 )V r 4(yU, r) - 6N c g 3 (fi)d(ji)5(r 



. L+ , L 



ldV 2 (n,r) 
c-{fjL,m 1 ,Tn 2 )-- 



r dr 

This is our new formula for the spin-dependent quark-antiquark potential. 

In Eq. (|33D, if we take each coefficient to be its tree level value, i.e., c 3 (/i,m) = 
c + (fi,mi,m 2 ) = 1 and c_(/i, mi,m 2 ) = d(fi) = 0, our result reduces to the EFG formula. 
Next we compare our result with that in the one-loop QCD calculation. First we see that 
V 5 (fj J ,mi,m 2 ) introduced in Ref. is not an independent function. If we take the one- 
loop values of c 3 (/z, m), c±(fi,m 1 ,m 2 ), and then calculate the correlation functions to 
one-loop, our formula (|33j ) reproduces all the logarithmic mass terms in Refs. @||. 

In conclusion, the leading logarithmic quark mass terms emerging from the loop contri- 
butions are explicitly extracted and summed up by matching the effective theory and the 
full theory and solving the renormalization group equation. The discrepancy appearing in 
the EFG results and one-loop calculation can then be understood. Our result shows that 
the effective theory can reproduce the full theory beyond tree level in 1/ m 2 and can be used 
in calculating the Green's functions with two heavy quark external lines. 
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